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THERMODYNAMIC  EXTRAPOLUT I ON  OF  ROCKET  ENGINE  PERFORMANCE 
PARAMETERS 

Ge  Mlnglong 

Abstract 

This  paper  uses  one  isentropic  reference  line  and  two 
isenthalpic  partial  derivatives  to  establish  extrapolation 
formulas  of  rocket  engine  performance  parameters.  These 
formulas  can  be  used  for  the  extrapolation  calculations  of 
specific  thrust,  characteristic  velocity,  nozzle-area  ratio 
and  the  thermodynamic  parameters  at  the  combustion  chamber  and 
the  nozzle  exit  when  there  are  changes  in  the  initial  enthalpy 
of  the  propellant,  combustion  chamber  pressure  and  nozzle 
exit  pressure. 

I.  Preface 

The  thermodynamic  calculation  method  for  rocket  engine 
ideal  performance  is  classic.  It  has  been  introduced  in  Refer¬ 
ence  1  and  other  books  and  periodicals  on  rocket  engines.  This 
diverse  calculation  is  very  useful  for  the  design  of  rocket 
engines  yet  the  time  expended  is  relatively  great  and  it  is  not 
convenient  to  arrange  and  put  into  book  form  the  large  amount 
of  data.  Therefore,  the  problem  of  the  extrapolation  calcula¬ 
tion  of  the  rocket  engine's  theoretical  performance  parameters 
was  brought  forth  early.  Many  people  have  given  formulas  for 
thermodynamic  extrapolation  calculation  [2-5]. 

lit  recent  years,,. following  ithe  popularization  of  pocket- 
electronic  calculators,  the  use  of  the  extrapolation  calcula¬ 
tion  method  to  determine  the  ideal  performance  of  rocket  engines 
is  even  more  convenient.  For  this  reason,  on  the  basis  of  ex¬ 
isting  extrapolation  formulas,  this  paper  derives  even  more 


accurate  extrapolation  formulas. 

II.  Isentropic  Reference  Lines 

Because  of  the  element  compositions  of  certain  propellants 
made  up  of  mixed  ratio  0/F  of  determined  value  are  fixed,  re¬ 
gardless  of  the  propellant's  initial  enthalpic  value  which  is 
affected  by  the  environmental  temperature  and  other  factors 

being  Ho  numerical  value  (H  or  H  )  and  no  matter  how  large 

co  c 

the  combustion  chamber  pressure  pc,or  pQ)  is  an<*  k°w  lar9® 

the  nozzle  exit  pressure  (p  ,p  , »p  „  p  ) ,  the  ideal  thermo- 

dynamic  process  of  the  rocket  engine  and  isentropic  limited 

balance  expansion  process  can  be  expressed  as  the  different 

isentropic  lines  c  t  e  e"...  etc.  on  the  same  enthalpy-entropy 

o  o  o 

chart  as  shown  in  Fig.  1. 


(3)**  • 


Fig.  1  Enthalpy-entropy  schematic  diagram. 

Key:  (1)  Enthalpy;  (2)  Isentropic  reference  line; 

(3)  Entropy. 

By  carrying  out  accurate  thermodynamic  calculations  of  a 
given  propellant's  initial  enthalpic  value  (HcQ)  and  combus¬ 
tion  chamber  pressure  (p_  ),  we  can  determine  the  thermodynamic 
parameters  of  combustion  chamber  state  point  cq  and  throat 

state  point  t  .  Afterwards,  we  selected  a  certain  noszle  exit 
o 


pressure  end  based  on  the  isentropic  conditions  we  calculated 

the  thermodynamic  parameters  of  these  nossle  exit  state  points 

(eo,e*,...).  On  the  isentropic  line  c0t0*0*" »  only  certain 

points  have  given  parametric  values  such  as  pressure  ratio  p  /p, 

c 

pressure  p,  temperature  T,  enthalpic  value  H,  molecular  weight 

M,  specific*  heat  at  constant  pressure  cp  etc.  and  the  values  of 

partial  derivatives/ <9  lnM^  and  /£lnM  \  (see  references  1 

( aTHp/  T  ITTH*  /p 

and  4  for  the  calculation  formulas  and  methods  for  these  two 
partial  derivatives).  Therefore,  this  isentropic  line  can  act 
as  the  reference  line  for  extrapolution  calculations.  A 
typical  isentropic  reference  line  is  listed  in  Table  1.  The 
parameters  on  the  top  half  of  the  table  were  taken  from 
Reference  5. 


m 


Table  1  Parameters  on  the  isentropic  reference  line  (gaseous 
hydrogen  and  gaseous  oxygen  propellant,  mixing  ratio 
O/F-7.9365,  initial  enthalpic  value  Hco*0  major 
calories/kg) . 

Key:  (1)  Parameter;  (2)  Conbustion  chamber;  (3)  Throat; 
(4)  Nozzle  exit;  (5)  Kg/cm2  (absolute;  (6)  Major 
calories/kg;  (7)  Kg/kg  mol;  (8)  Major  calories/kg; 

(9)  Data  obtained  by  extrapolation. 


On  the  isentropic  reference  line,  the  thermodynamic  para¬ 
meters  of  each  point  outside  the  known  parametric  points  can 
approximately  be  calculated  by  the  following  interpolation 
formulas : 


f-Om 

(AF),  -  £>(Aln,), 


Cl) 


In  the  formulas,  F  indicates  the  T,H,M. ..  etc.  parameters;  Q 
is  the  interpolation  coefficient  and  is  determined  based  on 
the  parameters  of  the  two  adjacent  known  points. 

For  the  isentropic  reference  line  in  Table  1,  the  calcula 
tion  results  of  QT,  QM  and  QH  etc.  interpolation  coefficients 
are  listed  in  the  lower  half  of  the  table. 


III.  Isenthalpic  Partial  Derivatives  and  Isenthalpic  Relational 
Formula 

Based  on  the  related  formulas  in  References  4  and  5,  the  two 
isenthalpic  partial  derivatives  used  in  this  paper  can  be 
indicated  by  the  following  formulas: 


In  the  formulas,  R  is  the  commonly  used  gas  constant  (1.98726 
major  calories/kg  mol  K) . 
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For  the  isentropic  reference  line  in  Table  1,  the  D^,  DM, 
and  values  calculated  by  each  of  the  above  formulas 

are  listed  in  the  lower  half  of  the  table. 


The  relational  formula  of  the  temperature  and  molecular 
weight  on  the  isenthalpic  line  is 


(AbT).  -  OrCAiv).  (4) 

UlalDa-DbCAfa*)*  (5) 


We  substituted  the  first  law  of  thermodynamics  expression 
and  the  state  equation  into  the  differential  formula  of  entropy 
and  obtained 


it 


JH  R  dp. 
T  U  9 


As  regards  isenthalpic  line  (dH=0) ,  from  this  formula  we  ob¬ 
tained 

(A*).- -£(4 Uf)m  (6) 

IV.  Establishment  of  the  Extrapolation  Formulas 


1.  Formulas  of  the  Combustion  Chamber's  Thermodynamic 
Parameters  and  the  Characteristic  Velocity 

When  the  propellant's  initial  enthalpic  value  Hc=Hco  and 

the  combustion  chamber  pressure  p_^P_^*  based  on  & H  =H  -H 

C  CO  c  c  o 

and  the  known  data  of  the  c  point,  we  can  calculate  the 

o 

parameters  of  the  c'  point  in  Fig.  1  from  formula  (1).  After¬ 
wards,  based  on  the  parameters  of  the  c'  point  and  p»p  ,  we 
can  calculate  the  parameters  of  the  c  point  from  formulas  (4) 
and  (5).  As  a  result,  we  obtained  the  combustion  chamber 
temperature  and  molecular  weight  formulas 
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T-  -  (r»  -  8^‘)  -*  !(°~ + £  - ^-)]  <?) 

Because  the  characteristic  velocity  c*  is  in  direct  ratio  to 

aJ  T  /M  ,  the  extrapolation  formula  of  the  characteristic  vel- 
c  c 

ocity  can  be  written  as 

*•  -  ciJTMJH.T,.  (9) 

2.  Formulas  of  the  Nozzle  Exit's  (p=pc)  Thermo  amir 
Parameters 

We  should  use  formula  (6)  between  points  c  and  c  and  points 
e  and  e ' : 


fcfe-M*fcfe  1 

t.  M«>  ft 

When  we  substitute  formula  (1)  into  this  formula,  we  obtain 

fc.fe._l n.jfe  + - (10) 

P~  P-  u„  + 

p* 

In  the  formula,  Meo  is  the  approximate  value  of  Me,  and  it  is  a 

known  molecular  weight  of  the  p*=p  point  close  to  p  ,  on  the 

©p  c 

isentropic  reference  line.  In  order  to  select  p  and  deter- 

eo 

mine  M^,  we  can  calculate  the  initial  value  of  pc,  according 
to  the  following  formula: 

*  ■  ... 

Key:  (1)  Initial. 

This  formula  was  approximately  obtained  from  formula  (10) .  The 

M  .  in  the  formula  is  the  known  molecular  weight  of  the 
© 


. .  v  4  ■ 
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isentropic  reference  line's  lowest  point. 


interpolation  coefficients  QT,  QDT,  QM,  QDM  and  0R  in  form¬ 
ulas  (7) ,  (8) ,  (10)  and  (11)  use  the  numerical  value  close  to 
PCQ.  When  ^Hc  <  0,  we  use  the  mean  value  of  c  -t  ;  when 
■A Hc  >  0,  we  use  the  numerical  value  of  p  ^  Pco* 


P_  i 

Based  on  the  known  data  of  the  e  point  and  the  value 

°  Pco 

calculated  from  formula  (10) ,  we  can  calculate  the  parameters  of 
e*  from. formula  (1).  Afterwards,  based  on  the  parameters  of 
the  e'  point  and  p=p ,  we  can  calculate  the  parameters  of  the  e 

C 

point  from  formulas  (4)  and  (5).  As  a  result,  we  obtain  the 


formulas  of  the  nozzle  exit's  temperature  and  molecular  weight 


T,  -  (T„  +  Prla^-)«p  [(  D*.  +  Q„bi&-)(hx  A-_b»  Ai)]  (12) 
y.  -  (*„  +  QJm&-)*p  [(D^  +  -  b  A-)]  (13) 

^4*  ^*4  ^44  ^4# 


3.  Formulas  of  the  Specific  Thrust  and  Other  Parameters 

■  When  we  substitute  H  =H  . =H  +  A H  and  H  =H  . =H  + 

Pc ,  c  c '  co  c  e  e '  eo 

QHln  — —  into  design  altitude  specific  thrust  formula 

*CO  r 

I«294 . 98  /y(Hc-He)/l, 000 ,  we  obtained 

I  -  9.32%  JH„  (14) 

interpolation  coefficients  Q^,  QDT,  Q^.,  QDM  and  QR  in  form¬ 
ulas  (12)  to  (14)  use  numerical  values  close  to  p  .  When 
PC'  .  co 

In -  \  O,  we  use  a  numerical  value  smaller  than  p  ;  when 

p„„  rco 

_  CO 

pc'  ^  O,  we  use  a  numerical  value  larger  than  PCQ* 

Pco 

The  nozzle  exit's  specific  area  and  the  nozzle's  area  ratio 
can  be  calculated  by  the  following  formula: 
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1.  -  U5iT./lfM.  (15) 

•  —  U*/‘*  (16) 

V.  Accuracy  of  the  Extrapolation  Formulas 

Based  on  the  above  formulas,  we  can  carry  out  extrapola¬ 
tion  calculations.  The  extrapolation  calculation  results  of 
the  hydrogen  and  oxygen  propellant  with  a  mixing  ratio  of 
0/F=7.9365,  nozzle  area  ratio,  characteristic  velocity,  combus¬ 
tion  chamber  temperature  and  designed  altitude  specific  thrust 
are  listed  in  Table  2,  3  and  4  and  Fib.  2.  In  order  to  make 
comparisons,  the  corresponding  precise  calculation  results 
[5,6]  and  the  extrapolation  calculation  results  of  Reference  5 
are  also  listed  in  the  tables  and  figure. 


*fp 
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000 
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w.-co  i  *+/* n 


2)aaa  6 (  ftwwii 

2.209  2.4C0  2.450  2.337 
7.130  7.10  7 .IN  7.270 
10.77  10.01  10.77  11.0 

14.70  14.01  14.70  15.01 

35.70  37.02  37.0  3«.45 
45.04  47.22  47.23  0.0 
0.49  M.H  C0.55  0.0 
70.0  0.21  04.07  0.0 
04.40  103.2  103.0  0.0 


Table  2  Comparison  of  nozzle  area  ratio  calculated  by  extra¬ 
polation  and  the  precise  values  (hydrogen  and  oxygen 
propellant,  mixing  ratio  0/F=7.9365,  combustion  chamber 
pressure  p  *6.805kg/cm2,  absolute). 

C 

Key:  (1)  Liquid  oxygen  and  liquid  hydrogen,  enthalpic 
value  Hc=-190.6  major  calories/kg;  (2)  Precise  value; 

(3)  Extrapolation  value  of  this  paper;  (4)  Extra¬ 
polation  value  of  Reference  5;  (5)  Gaseous  hydrogen 
and  gaseous  oxygen  enthalpic  value  Hc«629.1  major 
calories/kg;  (6)  Precise  value;  (7)  Extrapolation  value 
of  this  paper;  (8)  Extrapolation  value  of  Reference  5. 


9 


ftftftffft  » 

(1) 

/■**,  ft* 

IfflfWtt— * 

IfliM 

bS 

Mxmim 

HODCH 

□ 

*  ft  «(14) 

S3 

EaMffl 

I2EI 

5*X 

xms  (le 

00.35 

ar 

1 2197 

K22EZ3 

2150 

21S7 

2324 

2320 

2329 

40.13 

2104 

2104  |  2144 

2142 

2144 

2300 

2310 

2310 

0.005 

2UC 

1  21U  1 

2004 

El 

2251 

2254 

2271 

( 1 7  )*!  ***«**••»♦***/» 


Table  3  Comparison  of  characteristic  velocity  calculated  by 
extrapolation  and  the  precise  values  (hydrogen  and 
oxygen  propellant,  mixing  ratio  0/F=7.9365). 

2 

Key:  (1)  Combustion  chamber  pressure  pc»  kg/cm  , 
absolute;  (2)  Gaseous  hydrogen  and  gaseous  oxygen 
propellant,  enthalpic  value  Hc=0  major  calories/kg; 

(3)  Precise  value;  (4)  Extrapolation  value;  (5)  This 
paper;  (6)  Reference  5;  (7)  Liquid  hydrogen  and  liquid 
oxygen,  enthalpic  value  ^=-190.6  major  calories/kg; 
(8)  Precise  value;  (9)  Extrapolation  value;  (10)  This 
paper;  (11)  Reference  5;  (12)  Gaseous  hydrogen  and 
gaseous  oxygen,  enthalpic  value  Hc=629.1  major 
calories/kg;  (13)  Predate  value;  (14)  Extrapolation 
value;  (15)  This  paper;  (16)  Reference  5;  (17)  Note: 
in  the  table,  the  unit  of  the  characteristic  velocity 
is  meters/second. 


■ftSffft 

(1)  Fr/ft**»ft* 

1 .0B 

M4 

4.00 

aa 

».42 

40.03 

01.25 

01.50 

(2) 

an 

sir 

3217 
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(3)  ftXftftftft 
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(4)  ft: 


Table  4  Comparison  of  combustion  chamber  pressure  temperature 
calculated  by  extrapolation  and  the  precise  values 
(liquid  hydrogen  and  liquid  oxygen  propellant,  mixing 
ratio  0/F=7.9365). 

2 

Key:  (1)  Combustion  chamber  pressure  pc,  kg/cm  , 
absolute;  (2)  Precise  value;  (3)  Extrapolation  value 
of  this  paper;  (4)  Note:  in  this  table,  the  unit  of  the 
combustion  chamber  temperature  numerical  value  is  K. 


Fig.  2  Comparison  of  design  altitude  specific  thrust  calculated 

by  extrapolation  and  the  precise  values  (hydrogen  and 

oxygen  propellant,  mixing  ratio  0/F=7. 9365,  combustion 

chamber  pressure  p  =6.805  kg/cm2,  absolute). 

c 

Key:  (1)  Design  altitude  specific  thrust  I  seconds; 

(2)  Gaseous  hudrogen  and  gaseous  oxygen  (enthalpy 

value  Hc=629.1  major  calories/kg);  (3)  Liquid  hydrogen 

and  liquid  oxygen  (enthalpy  value  Hc=-190.6  major 

calories;kg) ;  (4)  Precise  value;  (5)  Extrapolution  value 

of  this  paper;  (6)  Extrapolation  value  of  Reference  5. 

(7)  Pressure  ratio  p  /p. 

c 

We  can  see  from  Tables  2  to  4  and  Fig.  2  that  the  precision 
of  the  extrapolation  formulas  in  this  paper  are  relatively  high 
and  can  be  suitably  used  for  the  extrapolation  calculations  of 
the  real  parametric  ranges  of  rocket  engines. 
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THERMODYNAMIC  CALCULATION  OF  THE  SYMBOLIC  FORMULA  FOR  ROCKET 
ENGINES 

Fang  Zhaokui 

Abstract 

A  symbolic  element  is  substituted  for  chemical  elements 
which  cause  practical  chemical  reaction.  With  the  use  of  the 
symbolic  formula,  a  new  calculation  method  has  been  set  up  to 
calculate  theoretical  performance  of  the  propellant.  It  will 
be  able  to  devise  a  general  computer  program  for  chemical 
equilibrium  of  rocket  engines  which  does  not  depend  on  prac¬ 
tical  chemical  elements. 

In  the  past,  China  widely  used  a  type  of  chemical  thermo¬ 
dynamic  calculation.  Its  major  feature  was  that  the  chemical 
reaction  formula  must  be  linked  with  the  concrete  chemical 
elements  of  the  propellant.  This  made  it  difficult  to  design 
general  programs.  The  thermodynamic  calculation  of  the  symbolic 
formula  proposed  in  this  paper  eliminates  the  above  mentioned 
difficulties. 

I.  Symbolic  Elements  and  Their  Molecular  Formula 

We  give  a  pair  of  positive  integers  i  and  j,  define  a  set  of 
symbolic  elements  A^  for  any  propellants  with  j  type  elements 
and  which  produce  i  type  combustion  products  as  well  as  their 
corresponding  symbolic  element  atomic  number  a^.  We  use  the 
ordered  permutation  of  symbolic  element  A^ : 

(<•  (0- 

It  indicates  the  i  type  combustion  product  molecular  formula 
produced  after  propellant  combustion.  Symbol  j$T  indicates  the 
ordered  permutation  as  shown  in  formula  (1) . 


For  example,  we  use  the  symbolic  elements  A1#  k2,  k^,  k 4 
to  replace  the  four  element  propellant  of  C',  H,  0,  N  and  let 
the  i  type  combustion  product  be  C  0.  Naturally,  based  on  the 
definition,  its  symbolic  element  molecular  formula  is: 

« 

°°  -  n  **  (2> 

The  connection  of  the  i  and  j  naming  sequence  with  the  product 
and  element  is  arbitrary,  yet  as  soon  as  we  gave  the  relation¬ 
ship  formula  (2)  solely  determined  the  CO. 

II.  Chemical  Equilibrium  Equation  of  the  Symbolic  Element 
Combustion  Product 


Definition:  for  any  molecular  type  combustion  product 
Ajaji»  we  let  its  dissociation  reaction  combine  with  symbolic 


element  A^ : 


n  z  ^  (*■  ■  u  2*-*o 


(3) 


It  acts  as  the  chemical  equilibrium  reaction  formula.  Because 
it  is  invested  with  the  symbolic  element  effective  partial 
meaning,  it  coincides  with  the  laws  of  conventional  multiplica¬ 
tion: 
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Key:  (1)  If;  (2)  If. 


-x‘H, 


<4> 


From  the  chemical  equilibrium  reaction  theory,  the  chemical 
equilibrium  equation  corresponding  to  formula  (3)  is: 

/(**  *  ***  *  0  -1^  -i)  C3) 

In  the  formula,  separately  indicate  the  partial 

pressure  of  symbolic  element  molecular  productliA^a^  and 


atomic  product  in  the  combustion  process  and  f^(T)  is  the 
equilibrium  constant  o£  this  dissociation  reaction  under  tem¬ 
perature  T.  If  there  is  basically  no  certain  element  existing 
in  this  molecule  (for  example,  basically  no  hydrogen  content 
in  the  CO  molecule  ) ,  then  P°  =1  because  of  the  conventional 
mathematical  law  and  we  must  J  be  able  to  satisfy  equation  (5) . 
We  used  the  logarithm  for  formula  (5)! 

-  w,ct) + 2  «) 

By  using  formula  (6) ,  we  can  cause  the  molecular  type  partial 
pressure  to  transform  into  atomic  type  partial  pressure, 
eliminate  the  molecular  type  variable  and  thus  greatly  reduce 
the  number  of  equations. 

III.  The  Symbolic  Element  Propellant  Matrix  and  Its  Normalized 
Process 

Definition:  we  assume  there  is  a  complex  propellant  composed 
from  the  n  type  pure  chemical  substance  mixture  and  molecular 
formulas  of  each  type  of  pure  chemical  substance.  The  enthalpy 
’'M  mixture  weight  percentage  are  separately 
i  XI  *.  (•  -  1>  2,  ••••)  and  then  the  characteristics  of 

this  propellant  can  be  completely  determined  by  the  following 
matrix: 


Formula  (7)  is  called  the  propellant  matrix.  During  performance 
calculations,  it  is  only  n-cessary  to  consider  the  theoretical 
performance  of  the  unit  mass.  For  this  reason 

Definition:  we  take  the  symbolic  element  molecular  formula 


of  a  unit  mole  propellant  with  a  molecular  weight  of  1  as  the 
standardized  equivalent  weight  formula  of  the  propellant.  The 
process  in  which  it  transforms  from  the  original  propellant 
to  the  symbolic  element  propellant  equivalent  formula  is  called 
the  symbolic  element  propellant  standardization  process. 


Theorem:  we  assume  there  is  complex  propellant  (7) ,  is 


the  atomic  weight  of  symbolic  element  A  ^ .  When  we  let  matrix 
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then  the  normalization  process  of  propellant  (7)  is  realized  by 
the  matrix  operation  of  the  following  formula 


A(j  +  l)  -  +  0  x  *00  - 


(9) 


In  the  formula,  A'(n,j+1)  is  the  transposed  matrix  of  matrix 
A(n,j+1)  obtained  after  deleting  the  last  column  in  matrix  (7). 
After  normalization  of  (9)  there  is  a  column  matrix.  The  lower 
symbol  "0”  of  column  matrix  element  a^Q  indicates  that  it  does 
not  belong  in  combustion  product  number  "i”.  From  (9)  we  can 
deduce 

i  £«?»-' 

This  is  the  inevitable  result  of  standardisation  processing. 

For  double  element  liquid  propellants,  if  the  combusting  agent 
of  corresponding  formulas  (7)  and  (8)  is  A(n,j+2),  B(n) ,  the 
oxidizing  agent  is  A(m,j+2),  B(m)  and  K  is  the  mixing  ratio. 


then  the  standardization  process  of  this  type  of  complex 
double  element  propellant  is: 

JQ  +  1)  -  i  +  1)  X  *(.)  +  KA'{m,  j  +  1)  X  B(m) J/(l  +  K)  (11) 

Finally,  we  obtained  a  column  matrix  in  the  same  way.  The 
proof  of  (9) ,  (10)  and  (11)  is  deleted. 

IV.  Basic  Equations  of  the  Thermodynamic  Calculation  of  the 
Symbolic  Formula  and  the  Program  Design 

Assuming  the  propellant  is  determined  by  positive  integers 
i  and  j  and  normalized  matrix  A(j+1) ,  then  the  combustion 
reaction  standard  formula  is: 

i  •  • 

u  n  a** — *  2>n  02) 

The  above  formula  shows  that  when  the  M  mol  propellant  com¬ 
busts,  it  separately  produces  combustion  product  ir  A^^ 
with  partial  pressure  P^,  yet  when  there  is  a  condensed  phase, 
P^  should  be  viewed  as  the  molar  number  of  the  condensed  phase 
product.  Now,  the  four  basic  equations  of  the  thermodynamic 
calculations  of  the  symbolic  formula  can  be  rewritten  as: 

1.  Conservation  of  mass  equation: 

#  03) 

2.  Conservation  of  energy  equation: 

(«*-*«***)  (1) 

OMHMM**)  C2) 

Key:  (1)  Isenthalpy-combustion  chamber  condition; 

(2)  Inentropy-nozzle  condition. 
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3.  Chemical  equilibrium  equation: 


In  formulas  (14)  and  (15),  (H°)i#  ^ST^i  and  *FT*i  are  sePar“ 
ately  the  enthalpy,  entropy  and  free  energy  when  the  i  type 
product  is  in  temperature  T. 

:  fJHsP,  (l)  (3) 

/'"to  ti  ohtwmb)(4)  • 

Key:  (1)  For  the  gas  phase;  (2)  Tor  the  condensed  phase; 

(3)  For  the  gaseous  phase;  (4)  For  the  condensed  phase. 

R  is  a  molar  common  suitable  constant.  The  differentiating 
formula  of  the  solidifying  phase  is  deleted. 

4 .  Dalton ' s  law : 

ik-Sft  <«urn»  (I<> 

Key:  (1)  Only  for  the  gaseous  phase. 

Solving  the  above  mentioned  equations  when  A(n,j+2), 

A(m, j+2) ,  K,  PQ,  Ho  and  are  known  forms  the  basic  proposi¬ 
tion  of  the  thermodynamic  calculations. 

In  the  program  design,  we  must  linearize  the  group  of  equa¬ 
tions,  use  (6)  to  reduce  the  number  of  equations  and  use  the 
partial  derivative  approximation  and  convergence  factor  to 
solve  the  problem  of  the  arbitrariness  of  the  initial  value  in 
a  large  range.  He  used  BCY  language  to  draw  up  the  commonly 
used  program  for  the  thermodynamic  calculation  of  the  symbolic 
formula  and  the  calculations  carried  out  on  a  109C  computer  for 
several  tens  of  propellants  were  all  successful.  The  above 
contents  were  limited  by  space  and  we  will  not  go  into  further 
detail. 
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